TILE COUNT IN THE INTERIOR OF REGULAR 2n-GONS 
DISSECTED BY DIAGONALS PARALLEL TO SIDES 

RICHARD J. MATHAR 



Abstract. The regular 2n-gon (square, hexagon, octagon,. . . ) is subdivided 
into smaller polygons (tiles) by the subset of diagonals which run parallel to 
any of the 2n sides. The manuscript reports on the number of tiles up to the 
78-gon. 



1. Summary 

Given the N = 2n sided regular polygon, its interior is dissected into non- 
overlapping regions (polygons, tiles) by the n(n — 2) diagonals parallel to any of 
the polygon's sides that connect edges of the polygon. The number of unrestricted 
diagonals is n(2n — 3). n{n — 1) of these are not considered here, so the number of 
regions F remains smaller than those generated by all of them [H El E] . 

The N sides define N/2 different non-oriented directions, each represented by 
a bundle of N/2 — 2 parallel diagonals; the subtraction of 2 indicates the polygon 
sides are not included in the count. The tile counts are summarized in Tables [1] 
and H 

N 4 6 8 10 12 14 16 18 20 22 24 26 28 

_n 2 3 4 5 6 7 8 9 10 11 12 13 14 

tiles F 1 6 25 50 145 224 497 630 1281 1606 2761 3302 5265" 
edges E 4 12 48 80 276 378 960 1062 2500 2860 5424 5980 10388 

Table 1 . Common results of the manual count of Figures IT1 [T51 

and of the C++ program of the appendix. 



By the discrete rotational symmetry, the polygon is invariant if rotated by mul- 
tiples of the angle 2ir/N. One can count the nonequivalent tiles in a segment in 
one particular ray direction, then multiply by N, and increase the count by 1 to 
include the central tile if n is even. 



2. Illustrations 

To verify the counts of Table [TJ a visual inspection follows for the cases up to 
the 28-gon. The tiles are enumerated in two different colors to demonstrate that 
each is in exactly one of the 2n replicas. Each of the numbers 1, 2, 3, ... , [F/N\ 
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N 68 70 72 74 76 78 

n 34 35 36 37 38 39 

F 204545 217210 258265 273282 321937 338208 
E 408068 420840 515376 530432 642580 656526 
Table 2. Results of the C++ program complementing Table [1] 



appears N times in a figure. From Figure [10] on, the enumeration covers only three 
rays to simplify the inspection of crossing patterns in the cluttered areas. 

The standard coordinate system puts the polygon corners on the unit circle at 
Cartesian coordinates (cos(7rfc/n), sin(7rfc/n)) for k = 0, . . . , N — 1. If these can be 
expressed as square roots [2] [3l |H [8j Q] [10] , the coordinates of the diagonal crossings 
are also of this form, since only determinental combinations of the coordinates of 
the polygon corners are involved to determine their positions. 




Figure 1. N = 6 sides: 6 tiles, one triangular tile replicated 6 
times. Corner coordinates are (±1,0), (±1/2, ±-\/3/2). 



TILE COUNT OF REGULAR 2n-GONS 3 




Figure 2. N = 8 sides: 25 tiles, one octagonal tile in the 
center and 3 triangular or kite-shaped tiles replicated 8 times 
outside the center. The 8 corner coordinates on the unit cir- 
cle are (±1,0), (0,±1), (±v / 2/2, ±V2/2). Internal crossings oc- 
cur at (±1/2,±(V2 - l)/2), (±(1 - \/2),0), (0,±(1 - \/2)), 
(±(1 - >/2/2),±(l- V2/2) and (±(1 - y/2)/2, ±1/2). 
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Figure 3. N = 10 sides: 50 tiles, comprising 5 trian- 
gular tiles replicated 10 times. The 10 corner coordinates 
on the unit circle are (±1,0), (±(1 + %/5)/4, ±^2(5 - ^5)/4), 
(±(V / 5 - l)/4, ±a/2(5 + V 5 )/ 4 )- Internal crossings are at 
(±(3 - y/5)/4, ± 72^5^75 /4), (±1/2, ±^/5 - 2 ^/5)), (±( 5 - 
75)/4,±(V5- 2)^2(5 + V 5)/4), (±(1 + ^/4^^/2(5 - V5)/4), 
(0,±(3- V5)V 2 (5- V5)/4), (±(1- V5)/2,0) and (0,0). 
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N=12, 1+12*12=145 




Figure 4. N = 12 sides: 145 tiles. One tile in the center and 12 
tiles replicated 12 times. 




Figure 5. N = 14 sides: 224 tiles, comprising 16 tiles replicated 
14 times. 




Figure 6. N — 16 sides: 497 tiles. One tile in the center and 31 
tiles replicated 16 times after rotation. 
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FIGURE 8. N = 20 sides: 1281 tiles. One tile in the center and 64 
tiles replicated 20 times. The lower graph shows a zoomed section 
of the upper graph. 
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Figure 9. N — 22: 1606 tiles, consisting of 73 tiles replicated 22 times. 




Figure 10. N = 24: 2761 tiles, one in the center and 115 tiles 
replicated 24 times after incremental rotation. 
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Figure 11. N = 26: 3302 tiles, 127 tiles per ray replicated N times. 




Figure 12. N = 28: 5256 tiles, one in the center and 188 repli- 
cated N times. 
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Appendix A. C++ Line Intersection Program 

Another way of counting is to enumerate all E edges and all V intersections 
(vertices) and to use the planar version of Euler's formula 

(1) F=l + E-V 

to count the faces of the graph. (The presence of the 1 indicates that the area 
outside the polygon is not counted.) The program starts from the edges of the 
exterior lines of the polygon that connect the roots of unity in the complex plane, 
plus the edges of the diagonals parallel to any of these edges. This base set is then 
rescanned producing a working set of chopped edges by looping over the segments 
of the base set, 

(1) each time computing all points of intersection (including "touches") with 
any element of the working set, 

(2) if this splits the element of the working set, replacing the element in the 
working set by its fragments, 

(3) adding the base set element or (if intersected or touched) its sub-components 
to the working set. 

The loop finished, E becomes the number of elements in the working set. V is 
computed by gathering all 2E terminal points of these in a point set, discarding 
duplicates that have a mutual distance smaller than some noise threshold set by 
the floating point arithmetic. 

The program is simple; no attempt is made to take advantage of the TV-fold 
rotational symmetry of the graph. 

#include <iostream> 
#include <cmath> 
#include <vector> 
#include <set> 
#include <algorithm> 

using namespace std ; 

/** A point with two Cartesian coordinates in the plane. 
*/ 

class Point { 
public : 

/** Abscissa doordinate . 

*/ 

double x ; 

/** Ordinate coordinate. 
*/ 

double y ; 

/** Default ctor. Point at infinity. 
*/ 

PointC) : x(INFINITY), y (INFINITY) 

{ 

> 

/** Ctor from two known Cartesian coordinates. 

* @param xcoo The x coordinate. 

* @param ycoo The y coordinate. 
*/ 

PointCconst double xcoo, const double ycoo) : x(xcoo), y(ycoo) 

{ 

} 

/** Distance to another point. 

* @param oth The other point to compare with. 

* ©return The Euclidean distance between this point and oth. 
*/ 
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double dist (const Point & oth) const 
{ 

return hypot ( x-oth.x, y-oth.y) ; 

} 

/** A distance somewhat larger than the estimated floating point error in the 

* computation of point positions. 

*/ 

static const double pointFuzzy = l.e-10 ; 
protected: 
private : 

> ; /* Point */ 



/** Comparison of two points. 

* Two points are considered equal if their mutual distance is (in the current coordinate units) 

* smaller than the constant value defined by Point :: pointFuzzy . 
*/ 

bool operator== (const Point &a, const Point b) 
i 

return ( a.dist(b) <= Point :: pointFuzzy ); 

> 

/** Comparison of two points. 

* This operator supports operations on point sets of the STL, and the underlying ordering 

* is not relevant for this program here. 
*/ 

bool operator< (const Point &a, const Point & b) 

-c 

if ( a.dist(b) <= Point: : pointFuzzy ) 

return false; 
else if ( a.x < b. x-Point :: pointFuzzy) 

return true ; 
else if ( a.x > b.x+Point :: pointFuzzy) 

return false; 

else 

return a.y < b.y ; 

> 

/** A straight line segment. 

* A connection between two points. 
*/ 

class Line { 
public : 

/** The two terminal points of the segment. 
*/ 

Point pts[2] ; 

/** Ctor given a start and end point. 

* @param strt One of the two terminal points. 

* Qparam fin The terminal points at the opposite end. 
*/ 

Line (const Point k strt, const Point & fin) 
{ 

pts[0] = strt ; 
pts[l] = fin ; 

} 

/** Copy ctor. 

* @param oth The other line that defines the new instance. 
*/ 

Line( const Line &oth) 
{ 

pts[0] = oth.pts[0] ; 
pts[l] = oth.pts[l] ; 

} 

/** Euclidean length. 

* Qreturn Pythagorean distance between the two terminal points. 
*/ 

double lent) const 
{ 



RICHARD J. MATHAR 



return pts [0] . dist Cpts [1] ) ; 

> 

/** Determine the point on the line defined by parameter t. 

* @param t The parameter along the line in units of the segment length. Values of and 1 recall the terminal points. 

* Sreturn The point at pts[0] + t* (pts [1] -pts [0] ) . 

* If t was outside the interval [0,1], this point is outside the line segment, but still coplanar. 
*/ 

Point atParam(const double t) const 
{ 

return PointC t*pts[l].x +(1 . -t) *pts [0] . x , t*pts [1] .y + ( 1 . -t) *pts [0] . y ) ; 

} 

/** Detect intersection with another line. 

* @param[in] oth The other line to intersect with. 

* @param [out] tOth The point of intersection parametrized by oth. strt+tOth* (oth . f in-oth . strt) . 

* ©return the point of intersection given by the parametrization t in strt+t*(f in-strt) . 

* If the lines are parallel, a value of infinity is returned. This line segment here intersects with 'oth' 

* if 0<=t<=l and 0<=t0th <=1 . Touching line segments have one of these two parameters in [0,1] and the other 

* one very close to or 1 . 
*/ 

double inters(const Line &oth, double k tOth) const 

{ 

/* Computation solves for the two linear equations for the point of intersection, and looks at the 

* determinant to detect parallelism. 

*/ 

double A [2] [2] ; 
double rhs[2] ; 

A[0] [0] = pts[l].x - pts[0].x ; 
A[0] [1] = oth.pts[0].x - oth.ptsEl] .x ; 
A[l] [0] = pts[l] .y - pts[0] .y ; 
A[l][l] = oth.pts[0].y - oth.pts[l] .y ; 
rhs[0] = oth.pts[0].x - pts[0].x ; 
rhs[l] = oth.pts[0].y - pts[0].y ; 

/* determinant of the coefficient matrix A of the two linear equations */ 
const double det = determinant (A [0] [0] ,A[0] [1] ,A[1] [0] ,A[1] [1] ) ; 

/* Use some measure of the determinant relative to the length 

* of the two vectors to put a threshold in a floating point environment. 
*/ 

if ( fabs( det) < pointFuzzy*f abs (lenO * oth.lenO) ) 
{ 

tOth = INFINITY ; 
return INFINITY ; 

} 

else 
{ 

/* If the determinant does not vanish, Cramer's rule computes the two line paremeters t and tOth. 
*/ 

tOth = determinant (A [0] [0] ,rhs[0] , A [1] [0] ,rhs[l]) /det ; 
return determinant (rhs [0] , A [0] [1] ,rhs [1] , A [1] [1] ) /det ; 

} 

} 



/** Compute the two line segments if this line is split at parameter t 

* Qparam t The line parameter in the range (0,1). 

* ©return The two line segments. The first one representing the (0,t) interval of points, 

* the second one representing (t,l). 
*/ 

vector<Line> splitAt (const double t) const 
{ 

vector<Line> out ; 

if ( inSplitRange(t) ) 

{ 

/* create and attach the two line segments */ 
out .push_back( Line(pts[0] , atParam(t) ) ); 
out .push_back( Line (atParam(t) ,pts [1] ) ) ; 

> 

else 

/* invalid specification. Return one component only, the segment itself. 



TILE COUNT OF REGULAR 2n-GONS 



17 



*/ 

out .push_back( Line(*this) ); 
return out ; 

} 

/** Decompose the line into non-overlapping segments defined by a set of line parameters of the new terminal points. 

* @param t The vector with numbers between and 1 that define the locations of the splitting location. 

* ©return The segmented version. The union of all these segments is the line itself. 

* ©warn The current implementation does not check that all t [] are in the range [0,1]. 
*/ 

vector<Line> splitAt( vector<double> t) const 

{ 

/* start and end point are always a start and end point in one of the new segments */ 
t.push_back(0.0) ; 
t .push_back(l .0) ; 

/* Sort the parameters numerically in ascending order 
*/ 

sort( t. begin O, t.endO) ; 

/* Eliminate points that are the same, so are neighbors in the sorted list. 

* Work with the local copy: move upwards in the list and delete the points later in the list. 
*/ 

forC int refidx = ; refidx < t.sizeQ-1 ; ) 
{ 

if ( fabsCt [refidx] -t [refidx+1] ) < pointFuzzy ) 
-C 

if C refidx == ) 

t.eraseC t.beginO + refidx+1) ; 

else 

t . erase (t . begin O+ref idx) ; 

> 

else 

refidx ++ ; 

> 

/* Duplicates (pairs, triples,...) are now removed and pairs of adjacent parameters in the list 

* define the new line segments to be returned. 
*/ 

vector<Line> out ; 

forC int refidx = ; refidx < t.sizeQ-1 ; refidx++) 
{ 

const Line seg(atParam( t[refidx]), atParamCt [refidx+1] ) ) ; 
out .push_back(seg) ; 

} 

return out ; 

} 

/** Decide whether the parameter t that runs from at the start 

* of the line segment up to 1 at the end is in the range, allowing 

* for some jitter in the floating point representation 

* Qparam t The line parameter. 

* ©return true if the line parameter is inside the range (0,1). 
*/ 

static bool inSplitRange ( const double t) 
{ 

return ( t > pointFuzzy kk t < 1 . 0-pointFuzzy ) ; 

} 

/** Decide whether the parameter t is close to one of the two values marking the terminal points. 

* @param t The line parameter. 

* ©return true if the line parameter is close to or 1 . 
*/ 

static bool inEndRange( const double t) 
{ 

return ( fabs(t) < pointFuzzy II fabs(t-l.O) < pointFuzzy ) ; 

} 

protected: 
private : 

static const double pointFuzzy = l.e-10 ; 

/** Determinant of the four values within a 2 by 2 matrix. 



18 



RICHARD J. MATHAR 



* Oreturn the value of aOO*all - a01*al0 . 
*/ 

static double determinant (const double aOO, const double aOl, const double alO, const double all) 
{ 

return a00*all-a01*al0 ; 

} 

} ; /* Line */ 



/** A set of line segments 
*/ 

class LineSet { 
public: 

vector<Line> lines ; 



/** Intersect all lines in the set and split them such that 

* all crossings are turned into terminal points of newly generated line elements 

* ©return The new line set with no intersections left. 
*/ 

LineSet splitAtlnters () const 
{ 

/* The updated list of lines to be returned. The outbound set. */ 
LineSet split ; 



/* Take the existing line elements one by one and intersect them 

* with all elements of the outbound set. 

*/ 

forC int lidx =0 ; lidx < sizeO ; lidx++) 
{ 

if ( split. sizeC) == ) 
split += lines [lidx] ; 

else 
{ 

/* The parameter list where lines [lidx] will be cut. 
*/ 

vector<double> tpoints ; 

/* The indices of elements in split [] that have been split by intersection with the current lidx . 
*/ 

set<int> rmPoints ; 

/* The subsegments of split [] that are created by cutting with the current lidx. 
*/ 

LineSet chaff ; 

/* loop over the line segments of the outbound set */ 
for(int slidx = ; slidx < split. sizeO ; slidx++) 
{ 

double t2 ; 

/* do both lines intersect or touch ? */ 

double t = lines [lidx] . inters (split . lines [slidx] , t2) ; 

if C isinfCt) ) 

/* no action if parallel, since in the current problem these segments are distinct */ 
else if C Line : : inSplitRange ( t ) ) 

{ 

if ( Line :: inSplitRange (t2) ) 

< 

/* They do intersect: mark position t for splitting, mark element slidx 

* for removal, and create its spliced segments for subsequent inclusion. 
*/ 

tpoints .push_back(t) ; 

chaff += split . lines [slidx] . splitAt Ct2) ; 
rmPoints . insert (slidx) ; 

> 

else if ( Line : : inEndRange (t2) ) 

/* line in the outbound set touches the line at lidx: keep line in the 

* outbound set and mark position t for subsequent splitting of lidx. 
*/ 

tpoints .push_back(t) ; 
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else if C Line : : inEndRange C t ) ) 
{ 

if ( Line : : inSplitRange (t2) ) 
{ 

/* line lidx touches the line in the outbound set. Same replacement 
* methodology for the outbound set as above, but keeping lidx intact. 
*/ 

chaff += split . lines [slidx] . splitAt Ct2) ; 
rmPoints . insert (slidx) ; 

/* Remove the outbound segments that have been split. Work in place which 

* requires working in reverse numerical order (from the maximum in rmPoints to the minimum) 

* to keep the indices intact. 

*/ 

while ( ! rmPoints . empty C) ) 
{ 

set<int> :: iterator rmidx = max_element (rmPoints .beginO , rmPoints . end() ) ; 
/* remove the full line in the outbound set */ 
split . lines . erase (split . lines . begin () + *rmidx) ; 

/* remove its index, so the next iteration of the loop will get the next smaller 
* index. */ 

rmPoints . erase (rmidx) ; 

> 

/* Insert the lines generated by splicing those removed above (list may be empty) */ 
split += chaff ; 

/* Add the line lidx to the set, or its components if split markers had been created above. 
*/ 

if( tpoints . empty () ) 

split += lines [lidx] ; 

else 

split += lines [lidx] . splitAt (tpoints) ; 

} 

> 

/* return the updated list of segments, which has no crossings (only contacts) left */ 
return split ; 

> 

/** Number of lines in the set. 

* ©return number of edges in the current list. 
*/ 

int sizeO const 
{ 

return lines. size () ; 

> 

/** Number of vertices in the set. 

* ©return The number of vertices. This counts all terminal points, the duplicates 

* at vertices shared by one or more edges only counted once. 
*/ 

int vertices () const 

{ 

set<Point> v; 

for(int lidx =0 ; lidx < lines. sizeO ; lidx++) 
{ 

/* insertion of the two terminal points of this line. 

* The Point :: operator== above eliminates duplicates at that time. 

*/ 

forCint i=0 ; i < 2 ; i++) 

v.insertC lines [lidx] .pts [i] ) ; 

} 

/* report the order of the set (count of element) finally residing in the point set 
*/ 

return v.sizeO ; 

> 



20 



RICHARD J. MATHAR 



/** Number of faces. 

* ©return the count of tiles as computed by Euler's formula. 
*/ 

int f aces C) const 

{ 

return 1+size () -vertices () ; 

> 

/** Add another line segment . 

* @param 1 The new additional line segment. 
*/ 

LineSet k operator += (const Line & 1) 

{ 

lines . push_back(l) ; 
return *this ; 

> 

/** Add another bundle of line segments. 

* @param 1 The list of additional line segments to be included. 

* Owarn The algorithm does not check for duplications. 
*/ 

LineSet k operator += (const LineSet & 1) 

{ 

lines . insert (lines . end() , 1 . lines . begin () , 1 . lines . end() ) ; 
return *this ; 

> 

/** Add another bundle of line segments. 

* @param 1 The list of additional line segments to be included. 

* ©warn The algorithm does not check for duplications. 
*/ 

LineSet k operator += (const vector<Line> k 1) 

{ 

lines . insert (lines . end() , 1 . begin () , 1 . end() ) ; 
return *this ; 

> 

protected: 
private : 

> ; /* LineSet */ 

/** Regular polygon with an even number of edges. 
*/ 

class Polyg : public LineSet { 
public : 

/** Ctor with count of edges 

* @param n half of the number of edges. 
*/ 

Polyg(int n) 
{ 

/* Construct the edges along the perimeter on the unit circle. 
*/ 

forCint e=0 ; e < 2*n; e++) 
{ 

Point a(cos( (double) e*M_PI/n) , sin ( (double) e*M_PI/n) ) ; 
Point b(cos((l.+e)*M_PI/n) , sin( (1 . +e) *M_PI/n) ) ; 
*this += Line(a,b) ; 

} 

/* Add the diagonals parallel to sides. 
*/ 

for(int e=0 ; e < n; e++) 
{ 

int e2 = e+1 ; 

for(int k=l ; k < n-1; k++) 

{ 

Point a (cos ((double) (e-k) *M_PI/n) , sin( (double) (e-k) *M_PI/n) ) ; 
Point b (cos ((double) (e2+k) *M_PI/n) ,sin((double) (e2+k) *M_PI/n) ) 
*this += Line(a,b) ; 

} 
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protected: 
private : 

} ; /* Polyg */ 

/** Main program. 

* @param argv There is one command line argument, the parameter n (half the number of edges) 
*/ 

int main (int argc , char *argv[]) 
{ 

int n = atoi Cargv [1] ) ; 

/* construct the polygon and the diagonals */ 
Polyg p(n) ; 

/* Split all pairs of edges of this graph at all crossings */ 
LineSet 1 = p . splitAtlnters () ; 

/* report the number of edges, vertices and faces present */ 

cout << l.sizeC) << 11 edges " << l.verticesC) « 11 vertices 11 « l.facesO << 11 tiles\n" ; 
return ; 

} 
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